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ABSTRACT 

If the Galactic rotation speed at the Solar circle is ~ 200 km s~^ or smaller, which is 
supported by several recent studies, the rotation curve of the Galaxy could be declining in the 
outermost region. Motivated by this, we investigate the effect of such declining rotation curve 
on the estimate of the MACHO mass and the fractional contribution of the MACHOs to the 
Galactic dark halo. Using Hernquist and Plummer halo models instead of the standard halo 
model, we find that the MACHO mass could be significantly smaller than that for the standard 
halo case. In particular, there exists a certain set of halo parameters for which the MACHO 
mass is O.IMq or less and at the same time the MACHO contribution to the total mass of the 
halo is almost 100 %. This result indicates that a halo which consists solely of brown dwarfs can 
be consistent with both of the observed microlensing properties and the constraints from the 
rotation curve, provided the outer rotation curve is indeed declining. 

Subject headings: dark matter — Galaxy : halo — Galaxy : structure — gravitational lensing 
— stars : low-mass, brown dwarfs 

1. Introduction 

The observations of gravitational microlensing towards the Large Magellanic Cloud (LMC) have 
revealed that a significant fraction of the Milky Way's halo consists of solar or sub-solar-mass compact 
objects, which are called massive compact halo objects (hereafter MACHOs; Alcock et al. 1993, 1996, 1997, 
Aubourg et al. 1993). So far at least six microlensing events have been detected toward the LMC, and 
recent analyses have shown that the MACHO mass is about 0.5Mq (Alcock et al. 1997). The fractional 
contribution of the MACHOs to the Galactic halo mass is also estimated to be / ~ 40%(Alcock et al. 1997), 
but may vary between 15 and 100% (e.g., Kan-ya, Nishi & Nakamura 1996). 

The mass estimate of ~ O.SA/q is consistent with red dwarf stars and old white dwarfs as baryonic 
candidates for MACHOs, but both halos have serious difficulties. A red dwarf halo cannot be consistent 
with the observations since red dwarfs do not appear in the star counts of Galactic halo fields (Graff & 
Freese 1996). If MACHOs are old white dwarfs, the initial mass function (IMF) of the progenitor stars must 
be strongly peaked around ~ 2Mq (Adams & Laughlin 1997, Chabrier, Segretain & Mera 1996), being 
completely different from the present day IMF (e.g. Scalo 1986). Moreover, the halos of distant galaxies 
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must be brightened by progenitor stars of the present day MACHO white dwarfs, which contradicts with the 
number count of distant galaxies (Chariot & Silk 1995). The overproduction of helium, carbon and nitrogen 
is also a major problem for a white dwarf halo (Chariot & Silk 1995; Gibson & Mould 1997). Furthermore, 
the age of the white dwarfs which are cool and faint enough to be consistent with the negative result of 
HDF star count (Flynn, Gould & Bahcall 1996) is very old (Chabricr, Segretain & Mora 1996). The age 
estimated by Chabricr et al. (1996) is not consistent with recent observations of the Hubble constant. 

On the other hand, brown dwarfs are possible and attractive candidates for MACHOs, since it is 
free from problems concerning metals and star counts. The largest problem of a brown dwarf halo is the 
consistency with the mass estimate of 0.5Mq. The mass estimate is, however, based on the 'standard' 
halo model, in which the rotation curve is flat out to the LMC. Hence, if the standard halo model does not 
represent the Galactic halo well, the MACHO mass would be changed considerably. Recently, Honma & 
Sofue (1996; 1997a) revealed that the outer rotation curve of the Galaxy could be declining if the Galactic 
constant ©o, which is the rotation velocity of the Galaxy at the Solar circle, is 200 km or less, and 
several studies indeed suggested O of ^ 180 km s"-'^ (e.g., Kuijken & Trcmaine 1994; Oiling & Merrifield 
1998). Furthermore, many observations have revealed declining rotation curves in the outer region of spiral 
galaxies (e.g., Casertano & van Gorkom 1991; Oiling 1996; Bland-Hawthorn et al 1997; see also Honma & 
Sofue 1997b and references therein), indicating that declining rotation curves are not uncommon. 

For these reasons, it is important to re-analyze the MACHO mass with 'non-standard' halo models for 

which the outer rotation curve decline. In this Letter we study the MACHO mass and the halo MACHO 
fraction based on such halo models and investigate the possibility of a brown dwarf halo. 



Since this investigation is motivated by declining rotation curves, we focus on the halo models whose 
density profiles decrease faster than the standard halo model, p oc r~^, of which Hcrnquist model and 
Plummer model are common examples. The density of Hernquist models is proportional to r"** at large 
radii, and it is known to trace the luminosity profile of elliptical galaxies moderately well. The density of 
Plummer models decreases faster than Hernquist model at large radii, proportional to r~^. The Plummer 
model is known to trace the light distribution of globular clusters to some degree. Conveniently enough, the 
potential-density pairs and the distribution functions are obtained analytically for both models. Below we 
summarize the analytic form of the potential-density pairs and distribution functions. 



2. 



Halo Models 



2.1. Hernquist Model 



The relative potential and the density for Hernquist model (Hernquist 1990) are written as. 
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where G is the gravitational constant, M is the total mass of the halo, and a is the scale length. Hernquist 
(1990) obtained the distribution function for this model as 
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where q = [aE/GMy/'^ and Vg = {GM/a)^/'^. Note that the relative energy £ is defined as f = * - 



2.2. Plummer Model 

The potential-density pair for the Plummer model (see Binney & Tremaine 1987) are 

GM 

(^2+^2)1/2' W 

and 

P^=[-^)ir'+^r'^'- (5) 
The distribution function for this model is written as 

Note that both of Hernquist and Plummer models contain only two parameters, the total mass M and 
the scale length a. 



3. Simple Estimate of r and m 

Based on the models described above, we calculate the optical depth r and the MACHO mass m. 
Paczynski (1986) defined the optical depth as 



m 

where Rb is the Einstein ring radius, m is the MACHO mass, and D is the distance along the line of sight. 
We assume that the distance to the LMC is 50 kpc, and also that all MACHOs have a unique mass m. r 
is related to the halo MACHO fraction, but independent of the MACHO mass because m cancels out in 
eq.(7). A quantity related to m is the differential event rate dV/dt (Griest 1991; see also Alcock et aL1996), 
which may be written as 



dV 327rttT f f R 



II 



, , FdVrdD. (8) 
dt m J J ^ ' 

Here wt is the maximum impact parameter, t = 2Re/v± is the Einstein ring crossing time, and Vr and 
v± arc the radial and tangential velocities of MACHOs, respectively. We set ut = 0.661 following Alcock 
et aL(1997). The total event rate F can be obtained by integrating eq.(8) with respect to t, namely, 
r = J (dT/dt) dt. The expectation of t can be obtained as 
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Note that < t > is proportional to y/m (Griest 1991), and hence, the MACHO mass can be estimated by 
comparing < i > with the observed one, < i >obs- To be conservative, we adopted Tobs — S.l^Qy x 10^^ 
and < i >obs of 80 days from six events so far detected (Alcock et al. 1997). Note that this < i >ohs is 
the simple arithmetic mean of the blending-corrected timescales for the six events, and the observational 
efficiency is not included here (see next section for more elaborate analyses). 

Figures 1 show the contours for r and m for Hernquist models and Plumnier models in the parameter 
space of a and M . The distance to the Galactic center i?o is assumed to be 7.5 kpc. If the lAU standard 
Rq of 8.5 kpc is assumed, the optical depth becomes slightly smaller, but the results are not significantly 
changed. Figures 1 demonstrate that the MACHO mass can be O.IM© or less depending on the halo 
parameters. In particular, the region for m < O.lA'/0 and the region for r ^ 2.1^0^ x 10~^ overlap 
(shadowed region in figures 1). This fact implies that the MACHO mass can be less than O.IMq while 
keeping the halo MACHO fraction close to ^ 100%. Therefore, a halo which consists solely of brown dwarfs 
(M < O.O8M0) is possible as long as these 'non-standard' halo models are accepted. 

Figure 2 shows the rotation curves modeled with the non-standard halos plus a disk and a bulge. An 
infinitely-thin exponential disk is assumed for the disk component, and Plummcr model is assumed for the 
bulge. The scale lengths are fixed to be 3 kpc for the disk, and 0.4 kpc for the bulge, respectively. The two 
rotation curves in figure 2 (full lines) correspond to the Galaxy's models with Hernquist and Plummer halo. 
For the Galaxy's model with Hernquist halo, we set M = 2.0 x lO^^M© and a = 13 kpc for the halo, and 
the disk mass of 3.5 x 1O-'^"M0 and the bulge mass of 0.8 x lO^^AfQ (hereafter model H). For the Galaxy's 
model with Plummer halo, we set M = 1.1 x IO^^Mq and a = 12 kpc for the halo, and the disk mass of 
4.0 X 10^°Mq and the bulge mass of 1.0 x IO^^Mq (hereafter model P). The halo parameters for models H 
and P lie in the shadowed region in figure 1. 

In figure 2 the observed rotation curve taken from Honma & Sofue (1997a) is also plotted (60 = 180 

km s""'^ assumed). The model rotation curves in figure 2 reproduce well the general tendency of observed 
rotation curves, except for a local deviation in the vicinity of the Sun. The local density of the exponential 
disks assumed here are 51 and 58 Mq pc^^, being in agreement with the studies on the local disk density 
(Kuijken & Gilmore 1989; Bahcall, Flynn & Gould 1992). Therefore, non-standard halo models which allow 
100 % brown dwarf halos does not conflict with the constraints from the rotation curve or the local disk 
density. This conclusion is not affected by small changes in the halo parameters. 

Why does the MACHO mass become so small when the non-standard halos are considered ? The 
substellar MACHO mass is obtained mainly because of the small velocity dispersion in the non-standard 
halos. The timescale t, which is the only observable used for the mass determination, is given by 
t = 2Re/v±. The Einstein ring diameter is proportional to \/?7i, and hence, m oc if the timescale is 
given observationally. Consequently, the MACHO mass becomes small if the velocity dispersion in the halo 
is sufficiently reduced. In fact, the tangential dispersion in the model P halo is 82 km s~^ at r = 20 kpc, 
which is less than half of the tangential dispersion for the standard halo, ^ 200 km s~^. 

4. Likelihood Analysis 

In this section we perform the maximum likelihood analysis and investigate more precisely the possible 
range of the MACHO mass and the halo MACHO fraction. The method applied here is basically the same 
as Alcock et aZ. (1997), whereas the halo models are the non-standard models described in section 2. The 
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likelihood of finding a set of A^obs detected events with duration times i — 1, Nohs{— 6), is written as ^ 

Nabs 



L(m,/) (iV(m))^°-=exp -fN{m)\ J| 
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where N{m) is the number of events expected from Poisson statistics (see eq.[8] in Alcock et al. 1997), / is 
the halo MACHO fraction, E = 1.82 x 10^ star -yr is the total exposure and e is the observational efficiency. 
For the observational efhciency we use the photometric efficiency of the MACHO year-1+2 observation 
given in figure 8 of Alcock et al. (1997). The likelihood contours for model H are shown in figure 3a, 
and those for model P in figure 3b. The estimated MACHO mass m is OMILomMq for model H, and 
Omto"otMQ for model P (the errors denote 68% confidence level). Most likely masses for both models are 
consistent with brown dwarfs. For the halo MACHO fraction /, the elongated contours in figure 3 make 
the possible range of / wide, but most likely values are fairly close to unity (0.92 for model H, and 0.86 for 
model P). The results for / as well as the MACHO mass m indicate that these halos can solely consist of 
brown dwarfs, confirming what we have found in the previous section. 



5. Discussion 

We have seen that a brown dwarf halo is a possible candidate for the Galactic halo when the 
non-standard halo models are considered. Assuming the halo MACHO fraction of unity for the models H 
and P, we obtain the MACHO mass within 50 kpc of 1.3 x IO^^Mq and 1.0 x IO^^Mq. These masses are 
somewhat lower than the estimate of 2lg ^ x IO^^A/q by Alcock et al. (1997), but do not contradict with 
each other when one considers the fact that the smaller Rq adopted here tends to make r higher for a 
given halo model. Since the total masses of model H and P galaxies are 2.4 x IQ^^Mq and 1.6 x 
70%~ 80% of the total mass may be in the form of brown dwarfs. To be consistent with the star counts 
at high Galactic latitude (e.g., Bahcall et al. 1994), such a brown dwarf halo should have a steep IMF at 
the low mass end (e.g., Kan-ya et al. 1996; Graff & Freese 1996). Unfortunately, any physical process for 
producing such a steep IMF is not known yet. However, the star formation from primordial gas may be 
very different from the present star formation because of the different cooling processes (e.g. Palla et al. 
1983, Uehara et al. 1996), and hence such an IMF may be possible if the halo brown dwarfs are different 
populations from disk stars. 

In order to test the possibility of a brown dwarf halo, direct observations of the halo brown dwarfs 
are definitely important. Unfortunately, star counts at high Galactic latitude so far made were not deep 
enough to detect the halo brown dwarfs. For instance, the star counts with HST (Bahcall et aZ.1994), which 
covers a sky area of 4.4 square arcmin with the limiting I-band magnitude of ^ 25 mag, detected no brown 
dwarf candidates. If the I-band absolute luminosity of m ~ O.O6M0 brown dwarfs is 18 mag (e.g., Delfosse 
et al. 1997), there should be ~ lO'^ brown dwarfs brighter than 25 mag in the vicinity of the Sun (model 
H assumed). The expected number in the HST filed of 4.4 square arcmin is, therefore, only 0.3, being 
consistent with the negative detection. For direct detection of the halo brown dwarfs, a deeper survey with 
a larger area is required. lOm-class telescopes with wide field of view, such as the SUBARU telescope with 
Sprime-Cam, will be useful for this kind of observation. 



The expression in eq. (10) is different from eq.{13) of Alcock et al. (1997) by the factor of (Ar(m))^°ba . Without this factor, 
L does not follow the Poisson distribution correctly. 
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Fig. 1. — Contours for the optical depth r (thin line) and for the MACHO mass m (dotted line) in the 

parameter space of M and a. (a) is for Hernquist model, and (b) for Plummer model. Contours for t 
correspond to 1.0, 1.5, 3.0x10"'' from left, and contours for m correspond to 0.04, 0.06, 0.14 Mq from 
left. In the shadowed area, t agrees with robs of 2.llg 7 x 10^'', and the MACHO mass is less than O.IMq. 



Fig. 2. — Model rotation curves for model H (thin line) and model P (thick line). Circles with error bars 
show the observed rotation velocities taken from Honma & Sofue (1997a). The Galactic constant 60 is 
assumed to be 180 km s""*^, and the radial distance is scaled with i?o=7.5 kpc. Rotation curves for each 
component are also shown. Dotted lines arc for disks, short dashed lines for halos, and long-dashed lines for 
bulges. See text for the model parameters. 



Fig. 3. — Likelihood contours of the MACHO mass m and the halo MACHO fraction / for model H (a) and 
for model P (b), respectively. The contours enclose total probabilities of 34%, 68%, 90%, 95%, and 99%. The 
crosses denote the peak of the likelihood. 
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